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ABSTRACT 

Recent studies have shown that during their coalescence, binary supermassive black holes (SMBHs) 
experience a gravitational recoil with velocities of 100 km s~ x uidck Si 600 km s — 1 . These velocities 
exceed the escape velocity v esc from typical dark matter (DM) halos at high-redshift [z <; 6), and 
therefore put constraints on scenarios in which early SMBHs grow at the centers of DM halos. Here we 
quantify these constraints for the most distant known SMBHs, with inferred masses in excess of 10 9 M Q , 
powering the bright quasars discovered in the Sloan Digital Sky Survey at z > 6. We assume that these 
SMBHs grew via a combination of accretion and mergers between pre-existing seed BHs in individual 
progenitor halos, and that mergers between progenitors with v esc < i>kick disrupt the BH growth process. 
Our results suggest that under these assumptions, the z ~ 6 SMBHs had a phase during which gained 
mass significantly more rapidly than under an Eddington-limited exponential growth rate. 
Subject headings: cosmology: theory - galaxies: formation - quasars: general - black hole physics 



1. INTRODUCTION 

The gravitational waves (GWs) emitted during the fi- 
nal stages of the coalescence of two merging black holes 
(BHs) carry linear momentum, implying that the cen- 
ter of mass of the system experiences a recoil (Bonnor 
& Rotenberg 1961; Peres 1962). The resulting recoil 
velocity had been known to be large, with values es- 
timated to be of order Wkick ~ 1000 km s^ 1 (Fitchett 
1983). Favata et al. (2004) has recently revisited this 
problem and computed recoil velocities, treating the spin 
and orbital dynamics of the merging BHs, as well as the 
generation of GWs in the strong gravity regime. They 
have found the range of possible recoil velocities to be 
100 km s _1 <J Wkick & 600 km s _1 , with the exact value 
depending on the mass ratio of the merging BHs, their 
spin, and orbital parameters. These velocities are large 
compared to the escape velocities of dwarf galaxies, and 
of the typical dark matter halos that existed at the early 
epochs of galaxy formation (z ^ 6). In hierarchical cos- 
mogonies, the SMBHs that are known to exist in the local 
universe grew via a combination of accretion and merg- 
ers between holes residing in individual DM halos. Mer- 
ritt et al. (2004) and Madau & Quataert (2004) recently 
considered several consequences of a large recoil that re- 
moves a BH from the center of its host galaxy. In partic- 
ular, Merritt et al. (2004) pointed out that the ejection of 
SMBHs from the shallow potentials of DM halos at high 
redshift implies a maximum redshift at which the progen- 
itors of present-day SMBHs could have started merging 
(and sticking) with each other. 

In this Letter, we consider the growth history of SMBHs 
that are in place at z > 6, and are thought to power 
the bright quasars recently discovered (Fan et al. 2000; 
2001; 2003) in the Sloan Digital Sky Survey (SDSS). As 



discussed in Haiman & Loeb (2001; hereafter HL01), rel- 
atively little time is available for the growth of these few 
xlO 9 M SMBHs prior to z ~ 6, and their seed BHs 
must be present as early as z ~ 10. A model in which 
stellar seed BHs appear in small progenitor DM halos is 
consistent with the presence of a ~ 4 x 10 9 M SMBH 
at z ~ 10, provided that each seed BH can grow at least 
at the Eddington-limited exponential rate, and that the 
progenitor halos can form seed BHs sufficiently eary on. 

The ejection of a merger-product SMBH from its host 
halo severely limits the ability of massive SMBHs to grow 
at z > 6, by disrupting the early stages of growth. In 
this Letter, as an example, we model the growth of the 
SMBH powering the most distant SDSS quasar, SDSS 
1054+1024 at redshift z = 6.43, with an inferred BH mass 
of ~ 4 x 10 9 Mq. Under the assumption that progenitor 
holes are ejected from DM halos with velocity dispersions 
c < ^kick/2, and do not contribute to the final BH mass, 
we find that typical recoil velocities must either be below 
the minimum value fkick = 100 km s _1 found by Favata 
et al., or else this SMBH must have had a phase during 
which it gained mass significantly more rapidly than the 
Eddington-limited exponential growth rate would imply. 

The rest of this Letter is organized as follows. In § 2, 
we discuss the inferred values of the relevant parameters 
(halo and BH mass) of SDSS 1054+1024. In § 3, we de- 
scribe the method we use model the growth of the SMBH 
by accretion and mergers. In § 4, we present our main re- 
sult, showing that excluding seed BHs from low-mass pro- 
genitor halos necessitates a faster-than Eddington growth 
rate. In § 5, we discuss various uncertainties about our re- 
sults. In § 6, we summarize the implications of this work 
and offer our conclusions. 1 



Throughout this paper, we adopt the background cosmological parameters for a flat universe as measured by the WMAP experiment, 
f2m = 0.29, Oa = 0.71, O;, = 0.047, h = 0.72, with a power spectrum normalization <Jg h -i = 0.9 and slope n = 0.99 (Spergel et al. 2003). 
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2. THE BLACK HOLE AND DARK MATTER OF HALO OF 
THE Z = 6.43 QUASAR SDSS 1054+1024 

The starting point for the constraints we derive below 
is simply the existence of a BH of mass Mbh at redshift 
z, residing within a dark halo of mass Mhaio- While the 
masses of SMBHs at the centers of nearby galaxies can be 
directly estimated, the mass of the SMBH powering the 
z = 6.43 quasar SDSS 1054+1024 is inferred indirectly 
from its observed luminosity. Under the assumption that 
the quasar emits a fraction rj of the Eddington luminosity, 
and using the template spectrum of Elvis et al. (1994) to 
make a bolometric correction, in the case of SDSS 1044- 
0125, we find M bh = 4.6 x lO 9 r?- 1 M . Given a sufficient 
fueling rate, bright quasars would naturally shine at their 
limiting luminosity, and we expect rj = 1. There is no ob- 
vious signs of beaming or lcnsing in the spectrum of this 
quasar (Willott et al. 2003). Indeed, the large observed 
size (~ 6 comoving Mpc; Mesinger & Haiman 2004) of 
its Strbmgrcn sphere makes it unlikely that the appar- 
ent flux of this quasar was significantly boosted by either 
lensing or beaming (Haiman & Cen 2002). In most con- 
ventional accretion models, r\ < 1, and the fiducial value 
of Mbh = 4.6 x 10 9 M Q would be a lower limit to the actual 
BH mass. However, there are models with 77 > 1; for ex- 
ample, in the recent "photon bubble" model of Begelman 
(2002), T) can be as high as r\ ~ 10, reducing the inferred 
BH mass. 

We next require the mass of the halo in which the quasar 
SDSS 1054+1024 resides. As in the case of the BH mass, 
the halo mass, or the velocity dispersion, can be directly 
measured for some nearby AGNs, but for distant quasars, 
we have to rely on indirect estimates. As described in 
HL01, Mh a io can be estimated based on the abundance 
of dark matter halos. In the case of the SDSS quasar, 
one bright quasar was found within a « 2000 deg 2 sur- 
vey area. In order to match this abundance, we find the 
halo mass has to be Mhaio ~ 8.5 x 10 12 M Q , with a corre- 
sponding velocity dispersion of 420 km s _1 . This result is 
weakly sensitive to the assumed duty cycle of quasar ac- 
tivity. We here assume £q = 4 x 10 7 (e/0.1)77 _1 yr, where 
e is the usual radiative efficiency of accretion. We use the 
standard Press-Schechter mass function (but find that our 
results would change little if we had instead adopted the 
recent numerical mass function in Jenkins et al. 2002, or 
the improved semi-analytical mass function of Sheth & 
Tormen 1999). Further details of the BH and halo mass 
determination are given in HL01. 

A relation between BH mass and halo circular velocity 
was recently determined in a sample of nearby galaxies 
(Ferrarese 2002). While this local relation does not neces- 
sarily hold at higher redshifts, it is interesting to note that 
it is in good agreement with the BH mass and (halo) veloc- 
ity dispersion our procedure yields for SDSS 1054+1024. 
Finally, the inferred halo mass is consistent with the value 
of Mhaio ~ 10 13 M derived from the spectral signatures 
of cosmic infall for this source (Barkana & Loeb 2003). 

3. GROWTH OF BLACK HOLES IN HIERARCHICAL MODEL 

In order to model the growth of the SMBH in our 
adopted ACDM cosmology, we rely on the merger his- 
tory of dark matter halos in the extended Press-Schechter 
(EPS) formalism (Press & Schechter 1974; Lacey & Cole 



1993). We follow HL01, and compute the central BH mass 
mirroring the assembly of its host halo. Given a parent 
halo of total mass Mhaio at redshift z, the EPS formalism 
specifies its average merger history back in redshift. Ev- 
ery branch of such a merger tree represents a progenitor 
of the parent halo, whose mass is continuously growing by 
accreting, and by merging with other progenitors. 

To keep our model simple, we assume that each pro- 
genitor of the parent halo develops a seed BH of mass 
Mgoed when the progenitor grows above a critical size, cor- 
responding to the velocity dispersion cr m i n = t>ki c k/2. Here 
Wkick is a typical recoil velocity for a coalescing SMBH bi- 
nary. Any seed BH that had appeared further up along 
the branch in the merger tree (e.g. the remnants of the 
first stars; Abel et al. 2002; Bromm et al. 2002) are thus 
assumed to be ejected from the progenitor halo, and not to 
contribute to the final SMBH mass at z ~ 6. This is a rea- 
sonable assumption, since in a typical "merger tree", each 
progenitor halo had been continuously undergoing merg- 
ers with other progenitors. The expected recoil velocity 
depends on the mass-ratio of the merging BHs; this de- 
pendence could be included in more detailed models of the 
merger tree, such as those utilizing Monte-Carlo realiza- 
tions (Kauffmann & Haehnelt 2000; Menou et al. 2001; 
Volonteri et al. 2003; Islam et al. 2003). Once the pro- 
genitor halo has merged with several others, its potential 
well will be deep enough so that the ejected seed BH could 
fall back into this enlarged halo. However, the ejected BH 
will likely have traveled far beyond the virial radius of 
the enlarged halo: we find tn u b(z)vki c k/ R v i r (z, a) <; 10 for 
z ~ 10 and a ~ 50 km s _1 . Madau & Quataert (2004) 
recently considered the dynamics of a recoiling BH in the 
fixed potential of the host galaxy, and showed that for 
a density distribution with a steep stellar cusp, dynami- 
cal friction can cause the BH to return to the halo center 
within ~ 10 6 yrs. It is unlikely, however, that steep stellar 
cusps exist in the earliest proto-halos that are thought to 
form a single massive star at z ~ 17 — 18 (which, as we will 
see below, is the birth redshift for the seeds responsible for 
the bulk of the final SMBH mass). 

We next assume that each individual seed BH 
subsequently grows exponentially by accretion, 
M hh (t) = exp[At(z, z/)/tacc]-Mseed, where i acc = 4 x 
10 7 (e/0.1)?7 _1 yr, as defined above, and At is the time 
elapsed between the formation of the seed BH at redshift 
Zf and a later redshift z. We assume that eventually, by 
redshift z, the smaller BHs in all progenitor halos coalesce 
together to form a single SMBH at the center of the par- 
ent halo (as long as the BH mergers are completed prior 
to redshift z, we do not need to specify when they take 
place). The mass of the resulting SMBH in the parent 
halo at redshift z is the sum of the individual BHs, each 
of which has grown by a different amount: 
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Here dP(z, z' , Mhaio, M) is the number of progenitors of 
mass M at redshift z' of a halo whose mass at redshift z is 
Mh a \o (Lacey & Cole 1993, eq. 2.15), and M m i n (z) is the 
mass of a halo whose velocity dispersion is Ukick/2- 

To summarize, our model for the assembly of BHs has 
five parameters. Two of these, Mhaio and 77, describe 
the observed quasar SDSS 1054+1024, and have relatively 
small uncertainties, as discussed in the previous section. 
The three parameters M soc( j, e, and Vkick relate to our 
model for the growth of the SMBH. The fiducial values 
of these parameters are chosen as follows. The seed mass 
is -Mscod = 10 M©, the typical value for a stellar remnant 
BH. VMOs (Carr, Bond & Arnett 1984) can leave larger 
remnants, weighing up to ~ 10 3 M© (Heger et al. 2003). 
The radiative efficiency is taken to be e = 0.1, based on 
the last stable orbit around a non-rotating BH. This value 
is consistent with a comparison of quasar light to remnant 
BH masses in nearby galaxies (Yu & Tremaine 2002; Allcr 
& Richstone 2002; Haiman, Ciotti & Ostriker 2004). A 
maximally rotating Kerr BH would produce a larger value 
of e = 0.42. 

For any given values of the above five parame- 
ters, equation 1 can be used to compute M b h = 
M bh (M" ha i , r;, M sccd ,e,w kick ). By requiring the predicted 
BH mass to equal the value inferred from observations, 
this relation can be inverted, and our model then yields a 
unique prediction for Vkick as a function of the five param- 
eters Mhaio, V, M sccd , e, and M bh . 



ingly higher rcdshifts z (shown explicitly in the bottom 
right panel) makes the contribution from the first ~ 20 
progenitors, forming at z <~ 18, dominant. 

No. of Seeded Progenitors Contribution to BH Mass 




4. RESULTS 

In our fiducial model, we find numerically that the max- 
imum recoil velocity that allows the growth of the SMBH 
in the quasar SDSS 1054+1024 is v kick = 64 km s" 1 . This 
value is significantly below the lowest values predicted by 
Favata et al. (2004) and Merritt et al. (2004). If actual 
recoil velocities are in excess of 100 km s -1 , this would be 
inconsistent with the fiducial SMBH growth model pre- 
sented here, and would require that some of the seeds grow 
their mass faster than the assumed Eddington rate. 

In order to illustrate the BH growth process in our model 
in somewhat more detail, Figure 1 shows the evolution 
of various quantities for SDSS 1054+1024 (this figure is 
an updated version of Figure 1 in HL01, which presented 
similar results for the earlier SDSS quasar 1054+1000 at 
z = 5.8). In this figure, we have assumed M scc( j = 10, e = 
0.1, rj = 1, and «kick = 64 km s _1 . With this combina- 
tion, equations 1 and 2 yield the required BH mass of 
M bh = 4.6 x 10 9 M Q at z = 6.43. The top left panel in Fig- 
ure 1 shows the number of progenitors of the parent halo 
(Mhaio ~ 8.5 x 10 12 M©) whose velocity dispersion exceed 
32 km s _1 . For reference, the bottom left panel shows 
the corresponding minimum halo mass. Going towards 
higher redshift, the number of progenitors increases, peaks 
at z « 11, and then decreases again as the typical progeni- 
tors are broken up into halos smaller than 32 km s _1 . The 
top right panel shows the contribution of progenitors from 
each redshift to the final black hole mass, and shows that 
the bulk of the BH mass is contributed by seed holes from 
17 <, z 18. There are no new seeds forming at z ^ 11, 
and the peak redshift is considerably higher than the peak 
at which most progenitors form. This is simply because 
the increased time available between z = 6.43 and increas- 



Fig. 1. — The assembly history of the black hole in the 
z = 6.43 quasar SDSS 1054+1024. The inferred BH mass 
is 4.6 x 10 9 M©, and the host halo mass is 8.5 x 10 12 M©. 
The four panels show, clockwise, the number of seeded 
progenitors (i.e. those with velocity dispersion above 
32 km s _1 ); the contribution of progenitors at different 
rcdshifts to the final BH mass at z = 6.43; the time 
available for the exponential growth of a seed between z 
and redshift of 6.43; and the halo mass corresponding to 
32 km s _1 at each redshift. 

5. DISCUSSION 

We have found above that, in order to grow an SMBH 
as massive as 4.6 x 10 9 M© in our fiducial model, we 
would need to utilize progenitors with velocity dispersions 
as small as a — 32 km s _1 . These small halos, how- 
ever, should be excluded from contributing to the final 
mass by the large recoil velocities. Another way of stat- 
ing our result is to note that in our fiducial model, but 
with er m i n = 50, 100, or 200 km s" 1 , corresponding to typ- 
ical recoil velocities of 100,200, or 400 km s -1 , the final 
SMBH mass at z = 6.43 is 5.2 x 10 s M©, 1.2 x 10 7 M©, 
and 1.2 x 10 5 M©, respectively - an order of magnitude or 
more below the inferred BH mass of SDSS 1054+1024. 

In order to assess the sensitivity of the result above 
to our assumptions, we here vary each of the parameters 
of our model. For each combination of Mh a i , r], M see di 
and e, we solve equation 1 with its left hand side set to 
Mbh = 4.6 x 10 9 ?7 _1 M©, and the halo mass set to Mh a i 
at z = 6.43, as discussed above. We further specify M see( j 
and e, and then find Wkick by a Newton-Rhapson method. 

We find that our results are insensitive to the adopted 
values of Mhaio and 77. Keeping all other parameters fixed 
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at their fiducial values, increasing or decreasing Mh a i by 
a factor of three yields i+ick = 82.2 km s _1 and ?+ick = 
49.6 km s -1 , respectively. This is not surprising, and re- 
flects the fact that the rare, massive halos at the tail of the 
mass function at z = 6.43 have similar merging histories. 
Similarly, increasing or decreasing r\ by a factor of three, 
we find Wkick = 80.4 km s _1 and fkick = 50.6 km s _1 , 
respectively. The sensitivity to the value of the final BH 
mass is only logarithmic because of the exponential growth 
predicted in equation 1. 

The sensitivity to the adopted value of the seed mass, 
M see( j, is the same as to 77 - changing the mass of the seeds 
or of the final BH is equivalent in our prescription, 2 since 
the model outlined above predicts the ratio M^/M scc d- 
However, the seed mass is more uncertain than the in- 
ferred BH mass. Once again keeping all the other parame- 
ters fixed at their fiducial values, we find that the choice of 
M S eed = 1, 10 3 , and 10 5 M Q results in w k i c k = 50.6, 157.8, 
and 332 km s _1 , respectively. Indications from recent 3D 
simulations (Abel, Bryan & Norman 2000, 2002; Bromm, 
Coppi & Larson 1999, 2002) are that the mass of the first, 
metal-poor stars are a few xl0 2 M Q . Nonrotating stars 
with masses between ~ 40 — 140 M Q and above <~ 260 M Q 
collapse directly into a BH without an explosion, whereas 
stars in the range <~ 140 — 260 M explode without leaving 
a remnant (Heger et al. 2003). Seeds as large as 10 5 M Q 
could only arise from the post-Newtonian instability of 
extremely massive stars. 

Finally, the largest sensitivity of our result is to the value 
of e. This is because e directly enters the e-folding time 
for the growth of the BH mass in equation 1. As an ex- 
ample, we find that the choice of e = 0.05, and 0.2 result 
in 438 km s _1 and 8.26 km s _1 , respectively. Thus, if the 
growth of each seed BH is as rapid as it would be with 
a typical radiative efficiency of ~ 0.05, the typical recoil 
velocities predicted by Favata et al. (2004) and Merritt et 
al. (2004) would still allow the build-up of the SMBH in 
the z = 6.43 quasar. Note that in this case, most of the 
final black holes mass at z = 6.43 would arise from the sin- 
gle most massive progenitor, starting to grow at z w 10. 
This is in contrast with our fiducial case with e = 0.1, 
which would imply that most of the BH mass was assem- 
bled by the addition of ~ 20 seeds, each of which started 
to grow at z <~ 18. A typical efficiency of ~ 0.05 would 
be significantly below the value obtained recently by Yu & 
Tremaine (2002). Yu & Tremaine (2002) compare the en- 
ergy density in quasar light and the mass density of local 
SMBHs (Soltan 1981) as a function of quasar luminosity. 
Interestingly, they find that bright quasars (with luminosi- 
ties similar to that of SDSS 1054+1024, or BH masses of 
J> 10 9 Mq) have a typical efficiency close to e ~ 0.2. This 
value is inferred only statistically, and therefore represents 
the average efficiency of the entire population of bright 
quasars at lower redshifts. Our results suggests that such 
a high value cannot hold for the past history of the indi- 
vidual source SDSS 1054+1025. 

The constraints obtained here could potentially be 
avoided if BH seeds at high redshift form only in a small 
fraction of the early dark matter halos (Madau & Quataert 
2004). In this case, when a dark matter halo, carrying 



a seed BH, merges with another dark halo, the seed BH 
would not experience a merger with another BH - it would 
then not be ejected, and it could continue to grow unim- 
peded. A small BH "occupation fraction" /bh would make 
the predicted abundance of high-z quasars rarer by a fac- 
tor of /bh, and the discovery of the z ~ 6 SDSS quasars 
could then be challenging to explain for /bh <C 1. How- 
ever, the fraction of DM halos harboring SMBHs would be 
expected to increase over time, and /bh could increase sig- 
nificantly between z <; 10 and z ~ 6 (Menou et al. 2001). 
More detailed models, following the merger history of BHs 
(Kauffmann & Haehnelt 2000; Menou et al. 2001; Volon- 
teri et al. 2003; Islam et al. 2003), are needed to quantify 
this scenario. 

Merging SMBHs at z <; 6 can loose a non-negligible 
fraction of their total mass to gravitational waves, espe- 
cially if they typically spin rapidly, and if they suffer a 
large number of mergers during their assembly history 
(Menou & Haiman 2004). More generally, although the 
mean expected recoil velocity is only weakly dependent on 
the BH spin (Favata et al. 2004), a high value for the 
typical spin of the merging BHs would strengthen our re- 
sults through (i) the loss of the BH mass to gravity waves, 
and (ii) through the larger expected value of the radiative 
efficiency e, increasing the accretion time-scale as e _1 . 

6. CONCLUSIONS 

The recent discovery of luminous quasars at redshift 
z > 6 provides evidence that supermassive black holes 
(SMBH) as large as several x 10 9 M Q were assembled dur- 
ing the first <, 10% of the current age of the universe. 
In the context of hierarchical structure formation scenar- 
ios, these early SMBHs grow via accretion and mergers of 
seeds that appear at much earlier epochs, z <; 10. Unless 
the growth by accretion of individual seed BHs is signifi- 
cantly faster than Eddington limited accretion at a fidu- 
cial radiative efficiency of e « 10%, we find that such sce- 
narios appear inconsistent with the large recoil velocities, 
100 km s _1 <, Wkick & 600 km s _1 , that were recently cal- 
culated to occur during the coalescence of massive BHs. 

A natural resolution of this discrepancy would be for 
individual SMBHs to grow in mass at significantly super- 
Eddington rates. In this case, the SMBHs can arise from 
seeds that appear relatively more recently, (i.e. closer to 
redshift z ~ 6, where the bright quasars exist), in DM ha- 
los whose potential wells are deep enough to retain these 
seeds despite their recoil. Forthcoming instruments, such 
as the James Webb Space Telescope (JWST) in the in- 
frared, or deep radio instruments, such as the Allen Tele- 
scope Array (ATA), Extended Very Large Array (EVLA), 
and the Square Kilometer Array (SKA), whose sensitivi- 
ties allow them to detect BHs as small as 10 5 M at z ~ 10 
(see Haiman & Quataert 2004 for a recent review) would 
be able to shed light on whether or not the seeds of z ~ 6 
quasar black holes indeed extend out to these redshift. 

The author thanks Kristen Menou and David Merritt for 
comments on an earlier draft of this paper, and gratefully 
acknowledges financial support from NSF grants AST-03- 
07200 and AST-03-07291. 



2 Note that the choice of r? and e have a small additional effect through the change they cause to the duty cycle, and hence to the inferred 
value of the halo mass. For sake of clarity, we chose here to neglect this effect, and parameterize the halo mass independently. 
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